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Terahertz (THz) radiation occupies a very large portion of the electromagnetic spectrum and has 
generated much recent interest due to its ability to penetrate deep into many organic materials 
without the damage associated with ionizing radiation such as x-rays. One path for generating 
copious amount of tunable narrow-band THz radiation is based on the Smith-Purcell free-electron 
laser (SPFEL) effect. In this Letter we propose a simple concept for a compact two-stage tunable 
SPFEL operating in the superradiant regime capable of radiating at the grating's fundamental 
bunching frequency. We demonstrate its capabilities and performances via computer simulation 
using the conformal finite-difference time-domain electromagnetic solver vorpal. 

PACS numbers: 29.27.-a, 41.85.-p, 41.75.Fr 
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Terahertz (THz) radiation is finding use in an in- 
creasingly wide variety of applications including medi- 
cal imaging, homeland security and global environment 
monitoring Increasing access to THz technologies 
requires the development of compact and tunable THz 
sources. Recent years have witnessed a resurgence of 
interest in Smith-Purcell free-electron lasers (SPFELs) 
operating as a backward wave oscillator ^ following 
on an idea initially discussed in Ref. [!]• The developed 
model was recently benchmarked by laboratory experi- 
ments [1, Q . THz sources based on SPFELs are foreseen 
to have table-top footprint and can operate in a contin- 
uous wave mode, enabling the production of moderate 
average output power (on the order of Watts) . 

In an SPFEL, a low energy (-- 50 keV) sheet DC 
electron beam is propagated close to a metallic grating 
with velocity v = cfiy . The beam excites evanescent 
surface waves with axial field of the form Ey,e{x,y) = 
£^o,e exp(aa;) where a = 2'K/{[3^\e) and 7 = (1 - /?^)^^/^ 
is the Lorentz factor. The evanescent wave can, un- 
der certain circumstances, have a negative group veloc- 
ity 0- such a case the wave counter-streams the 
electron beam direction and can couple to the beam, 
thereby giving rise to an energy modulation. Due to 
the non-relativistic nature of the beam (7 ~ 1), the 
impressed energy modulation eventually transforms into 
a density modulation at wavelength Ag. The produced 
microbunches will result in strongly enhanced radiation 
at harmonic frequencies of the microbunching frequency 
/e = c/Ag. In an SPFEL the radiative mechanism is the 
Smith-Purcell (SP) effect 0. If instead of a DC elec- 
tron beam, a beam consisting of a train of microbunches 
is used the SP radiation is emitted in the super-radiant 
regime (sl-fioj in which the radiation rate goes as th num- 
ber of electrons in each microbunch squared. Prebunch- 
ing the electron beam in a way that satisfies emission 
of super-radiant radiation is however challenging and a 
possible solution discussed in, e.g., Ref. [ll|] significantly 



decreases the average power capability of the SPFEL. 

In this Letter we consider and present detailed nu- 
merical simulations of a "two-stage" SPFEL; see Fig. [T] 
In such a configuration a first stage, referred to as a 
"buncher" , is optimized to enhance the beam-evanescent 
wave interaction, thereby resulting in faster bunch- 
ing than in the configuration analyzed in previous pa- 
pers A second stage, referred to as a "radiator" , 
consists of a grating with parameters tuned to produce 
coherent SP radiation at frequencies n/e {n is an integer). 




bunching stage radiating stage 



FIG. 1: (Color online) Diagram of the two-stage SPFEL. The 
rectangular box represents the 3D computational domain used 
in the Vorpal simulations. Grey and blue blocks respectively 
stand for perfectly matched layer and metallic boundaries. 

The numerical simulations were performed using VOR- 
PAL, a conformal finite-difference time-domain (CFDTD) 
particle- in-cell electromagnetic solver ^is']. vorpal is a 
parallel, object-oriented framework for three dimensional 
relativistic electrostatic and electromagnetic plasma sim- 
ulation. We extensively benchmarked our initial vor- 
pal simulations of the SPFEL process against earlier 
work [l3|; see Ref. [l^. The geometric configuration of 
our simulation model is shown in Fig. [T] The model in- 
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TABLE L Grating and beam parameters used for the Vorpal 
simulations. 



Parameter 


Symbol 


Value 


Unit 


Grating Period 




200 


lira 


Groove Width 


W 


100 


/im 


Groove Depth 


h 


100 


^m 


Number of Periods/grating 


No 


75 




Electron Energy 


E 


25-200 


keV 


Beam Current 


I 


> 135 A/m 


Beam Thickness 


b 


20-400 


^m 


Beam Clearance 


9 


20-30 


/im 


Grating gap 


2X0 


6 + 25 


^m 


External Magnetic Field 


Ba 


2.0 


T 



eludes perfectly-conducting rectangular gratings, a par- 
ticle source on the left end side and a "particle sink" on 
the right that allows macroparticles to exit the computa- 
tional domain without being scattered or creating other 
source of radiation. The grating grooves are along the z 
direction and the beam propagates along the y axis. Be- 
side the grating, the other surrounding boundaries con- 
sist of perfectly matched layers that significantly sup- 
presses artificial refiections of incident radiation. The 
particle source produces a uniformly-distributed (in all 
directions) DC beam with an instantaneous risetime and 
is transversely confined by an uniform external axial mag- 
netic field Ba = Bay. The parameters used for the sim- 
ulations presented below are gathered in Tab. HI with 
£: = 50 keV, / = 135 A.m-\ 6 = 50 /im, and 5 = 20 /im; 
for the sake of simplicity the buncher and radiator grat- 
ings are identical. 



Considering the two-dimensional problem, the evanes- 
cent waves produced from the beam passing over each 
grating located at ±a;o are Ey ,,{x, y) — £'o,e(y)e^"^'^"+'^' 
(where ± respectively designate the electric field contri- 
butions from the upper and lower grating). The over- 
lap of the evanescent waves from both gratings results in 
an evanescent field Ey^e{x,y) = 2EQ,e{y)e~"^" cosh{ax). 
Beside being stronger by a factor of two (for a; = 0) 
thereby bunching the beam stronger, the field is also 
more uniform across the beam's thickness and effectively 
result in a more uniform bunching. The strengthened 
evanescent field may also be taken advantage of to relax 
the start current requirements. These features are con- 
firmed by numerical simulation; see Fig. [2l The velocity 
modulation is approximately twice as large for the dou- 
ble grating configuration than for a single grating con- 
figuration. The presence of the second grating also af- 
fects the dispersion relation and Ae, which in-turn leads 
to a different bunching frequency than the single grat- 
ing system; see analysis in Ref. [16[. This feature offers 
a greater flexibility for tuning the evanescent wave fre- 
quency (and therefore the radiation frequency) by either 
varying the electron beam energy, or altering the gap be- 
tween the two gratings as demonstrated in Fig. |3l As 
the gap increases, the frequency converges with that of 
a single-grating sytem. In addition the double grating 
configuration allows the emission of SP radiation at the 
bunching frequency. Over the considered range of en- 
ergy E £ [25,200] keV and grating gap 2xo e [90,440], 
significant velocity modulation is observed and the corre- 
sponding evanescent wavelength variation of respectively 
15 % and ~ 35 % are approximately two orders of 
magnitude than the produced radiation spectrum width 
as will be shown later. 
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FIG. 2: (Color online) Snapshot of the axial phase spaces 
{y,Sy) (where Sy = jPy/lPyiy = 0)7(1/ = 0)] - 1 is the frac- 
tional axial momentum spread) for a double (a) and single (b) 
grating taken at 151 ps. The lower insets on each of the plot 
present zoomed-in phase space over for y G [0.50, 0.53] cm, for 
which the macroparticles have been color-coded accordingly 
to their height within the beam. The double-grating simula- 
tion has less distinct layers, signifying more uniform bunching 
across the height of the electron beam. The beam parameters 
are £; = 50 keV and I = 135 A.m'^ 

We first analyze the buncher section of the proposed 
two-stage SPFEL. Contrary to the SPFEL configurations 
explored in previous work, the beam propagates between 
two identical gratings arranged symmetrically; see Fig.[T] 



FIG. 3: (Color online) Evanescent wave frequency versus 
beam energy for the single-grating configuration (a) and 
evanescent wave frequency versus grating separation for a 
double-grating configuration (b). The markers are results 
from VORPAL simulations while the dashed line is obtained 
from the dispersion relation derived in Ref. 9]. In both 
graphs, the green solid line denotes the minimum allowed SP 
frequency for a grating period of 200 /im. 

The second stage of the proposed SPFEL consists of 
passing the microbunched beam close to a single grating. 
The spectral intensity radiated by a bunch of N electrons 
via the SP effect is related to the single electron intensity 
£li via £U =jh\iW + N^\S{u:)\^] where Siu;) is the 
intensity-normalized Fourier transform of the normalized 
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charge distribution S{t) [T^l- Considering a series of A'^f, 
identical niicrobunches with normahzed distribution A{t) 



we have S{t) = N^^ J2nU + nT) (where T = \Jc is 
the period) giving 15(^)1^ = ^|A(a;)p. The intra-bunch 
coherence factor ^ = N{;^ sm\ujNbT /2) /[sii? [ujT 12)] 
describes the enhancement of radiation emission at reso- 
nance, i.e. for frequencies u> = 27rnc/Ae. 
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FIG. 4: (Color online) Snapshots &t t = 906 ps of the axial 
electric field Ey and corresponding spatial distribution of the 
beam in the {x,y,z = 0) plane for a two-stage SPFEL (the 
colorbar units are MV.cm"^). The associated axial projec- 
tions at two axial locations are shown in (c) and (d). 

Figure |4] shows a contour plot of the E'y-field over 
the computational domain and clearly demonstrates the 
emission of SP radiation with a planar wavefront from the 
second grating. The planar wavefront is a signature of 
superradiance due to constructive interference between 
successive bunches. The wavevector associated to this 
radiation makes a 106° angle and the wavelength value 
~ 540 jiva. is consistent with the one exected from the SP 
relation Q Am — Ag/|m|(/3^-'^ — cos6') = 540 /im (tak- 
ing TO = 1, Ag = 200 //m and (3 = 0.4127). Whereas 
the single-grating SPFEL requires that the super-radiant 
emission occurs at m > 2, the proposed two-stage SPFEL 
can operate in the super-radiant mode at m = 1 by a 
proper selection of parameters. Even if the radiator and 
the buncher have the same grating period, the increased 
bunching frequency from the double-grating dispersion 
relation may match that of the m=l mode of the SP re- 
lation. Operating at to = 1 enables the generation of 
higher radiation rates compared to to > 2. Figure 0] 



also displays a snapshot of the spatial distribution and 
associated projections confirming the formation of mi- 
crobunches. 

The time evolution of the magnetic field recorded 
aX 6 = 106° and R = 2 cm along with its spectrum are 
presented in Fig. [5] The spectrum shows that signifi- 
cant emission occurs at fe = 556 GHz, the frequency 
of the evanescent wave supported by the buncher sec- 
tion, and its harmonics. The bandwidth of the radiation 
emitted at the first harmonic is Sf ~ 1.3 GHz result- 
ing in df/fe — 0.23 %. The spectrum also displays a 
413 GHz frequency component which corresponds to the 
evanescent wave produced by the downstream grating. 
Taking the wavefront to be a plane wave with steady- 
state peak magnetic field of Bq ~ 40 ^T (see Fig. [5]), 
we estimate the time-averaged Poynting vector to be 
(5) = cB^/{2fio) ~ 2 X 10* W.m-2. Harvesting the 
radiation using a 1 x 1 cm^ mirror would result in a total 
average power of P ~ 2 W. 
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FIG. 5: (Color online) Time evolution of the magnetic field 
at {R,6) = {2 cm, 106°) (a) and associated FFT (b). The 
beam parameters are I = 135 A.m~^ and E = 50 keV. The 
inset in (b) is a zoom- in of the 1st harmonic peak showing the 
bandwidth of 1.3 GHz (FWHM). The numbers (1), (2) and 
(3) indicate the harmonics of the evanescent wave associated 
to the double grating configuration. The green vertical arrow 
indicates the evanescent wave supported by the downstream 
grating, which is below the minimum frequency allowed by SP 
radiation and thus only scatters off the edges of the grating. 



In summary we have demonstrated several advantages 
of a two-stage SPFEL and show that the device can 
be used to produce tunable super-radiant THz radiation 
with average power of the order of Watts. 

We thank Dr. P. Spentzouris (FNAL) for granting 
us access to the NERSC computers where our simu- 
lations were performed. This work was partially sup- 
ported by the US Department of Education under con- 
tract P116Z050G86 with Northern Illinois University. 



[1] M. Tonouchi, Nature Photonics, 1, 97-105 (2007). 

[2] H. L. Andrews and C. A. Brau, Phys. Rev. ST Accel. 

Beams, 7, 070701 (2004). 
[3] V. Kumar and K.-J. Kim, Phys. Rev. E 73, 026501 

(2006). 



[4] J. M. Wachtel, J. Appl. Phys., 50 (1) (1979). 

[5] H. L. Andrews and C. A. Brau, Phys. Rev. ST Accel. 

Beams 12, 080703 (2009). 
[6] J. T. Donohue, L. Courtois, P. Modin, and J. Gardelle, 

Phys. Rev. ST Accel. Beams 12, 110701 (2009) 



4 



[7] S. J. Smith and E. M. Purcell, Phys. Rev. Lett. 92, 1069 

(1953). 

[8] H. L. Andrews, C.H. Boulwarc, C.A. Brau, J.D.. Jarvis, 

Proceedings of the 2004 PEL conference, 278 (2004). 
[9] D. Li, K. Imasaki, Z. Yang, and G.S. Park, Appl. Phys. 
Lett. 88, 201501 (2006). 
[10] A. Kesar, R. Marsh, R. Temkin, Phys. Rev. ST Accel. 

Beams 9, 022801 (2006). 
[11] Y. Li andK.-J. Kim, Appl. Phys. Lett. 92, 014101 (2008). 
[12] Z. Shi, Z. Yang, F. Lan, X. Gao, Z. Liang, and D. Li, 
Nuc. Instr. Meth. A 607, 367-371 (2009). 



[13] C. Nieter and J. R. Gary, J. Camp. Phys., 196, 448 
(2004). 

[14] J. T. Donohuc and J. Gardelle, Phys. Rev. ST Accel. 
Beams 8, 060702 (2005). 

[15] G. R. Prokop, P. Piot, M. G. Lin, and P. Stoltz, Proceed- 
ings of the 2009 PEL conference, 180 (2009). 

[16] H. P. Preund and T. M. Abu-Elfadl, IEEE Trans. Plasma 
Sci., 32, 1015 (2004). 

[17] J. S. Nodvick and D. S. Saxon, Phys. Rev. 96, 180 (1954). 



